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We study Thomson scattering from the antiferroquadrupole ordering phase in TmTe. On 
the basis of the group theoretical treatment, we classify the selection rules of the scattering 
intensity governed by the orientation of the scattering vector G. Then, numerical verification 
is performed by invoking the ground states which are deduced from a J = | multiplet model. 
The obtained intensity varies drastically depending on the magnitude and direction of G. 
We also calculate the scattering intensities under the applied field for H || (001) and (110). 
Their results behave differently when the orientation of G is changed, which is ascribed to the 
difference of their primary order parameters; 0% and 0\ for H || (001) and (110), respectively. 
We make critical comparisons between our results for TmTe and the experimental ones for 
CeBg. First, we assert that the intensities expected from TmTe at several forbidden Bragg 
spots are sufficient enough to be experimentally detected. Second, their intensities at (\\\) 
differ significantly and may be attributed to the difference of the order parameters between 
the IVtype (0\ and 0°) an d IVtype (O yz ,O zx , and O xy ) components, respectively. 

KEYWORDS: Thomson scattering, non-resonant X-ray scattering, multipole, antiferro- 
quadrupole order, forbidden Bragg spot, TmTe 

1. Introduction 

The interplay of orbital and spin degrees of freedom in localized magnetic materials brings 
about a wide variety of interesting phenomena. In many /-electron systems, due to the strong 
coupling between the spin and orbital angular momenta, the states are described by the 
multiplets of the total angular momentum J. When the symmetry exhibited by the system is 
sufficiently high, the multiplet enables the higher rank multipoles as well as the dipole (rank 
one) be active. In fact, various experimental and theoretical studies have been devoted to 
clarify the nature of the ordered phase of multipole order parameters with rank higher than 
one. 1,2 Among the most investigated systems, the materialization of the antiferroquadrupole 
(AFQ) ordering phase has been established in the materials such as CeE>6 and DyE>2C2- 

Among many experimental probes, scattering experiments such as resonant X-ray scat- 
tering (RXS) and (non-resonant X-ray) Thomson scattering provide very powerful tools to 
reveal the natures of the higher rank multipolar ordering phase. For example, the AFQ order- 
ing phases in CeB@ and DyE>2C2 are investigated in detail by means of RXS 3-7 and Thomson 
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scattering. 4 ' 8-10 

On the other hand, the situation of TmTe may still be rudimentary. This material is 
believed to show the AFQ order below Tq = 1.8 K. 11 Although many evidences for the AFQ 
order were gathered in terms of various experimental probes, 12-14 there remain some important 
issues unsettled yet. For instance, the crystal electric field (CEF) level scheme, the component 
of the primary order parameter, the nature of the multipolar interaction, and so on. In order 
to address such issues, the approaches in terms of the scattering experiments may be helpful. 
In our previous work, we have determined the CEF level scheme as T$ — Tq — T7 and analyzed 
some properties expected from the azimuthal angle dependence of the RXS intensity, which 
are useful to distinguish the type of the order parameter. 15 

In this paper, we carry out some investigations on Thomson scattering expected from the 
AFQ phase in TmTe. After introducing the theoretical framework to calculate the scattering 
intensity, we first proceed to classify the selection rules of intensity governed by the direction 
of the scattering vector. In the absence of the applied field, such selection rules as well as 
the domain consideration determine the whole intensity. The intensity exhibits the strong 
dependence on the magnitude and orientation of the scattering vector. We verify the quali- 
tative results with the numerical calculation performed on the theoretical model developed 
in our previous paper. 15 We also investigate how the application of the external field alters 
the scattering intensity. When the field is applied along (001) and (110), the primary order 
parameters derived from the model are 0° = 1/2(2 - J x - J y ) and 0\ = V3/2(J X — J y ), 
respectively. We find their intensities show different behaviors as a function of the orientation 
of the scattering vector, reflecting the difference of their primary order parameters. 

We also try to compare the present results with those obtained for CeBg. 4 ' 16 Although both 
TmTe and CeBg exhibit the AFQ ordering phases, it is said the components of the order pa- 
rameters are different from each other; the T3-type {0\ and O®) in the former and the Ts-type 
(O yz , O zx , and O xy ) in the latter. Here, O yz = V3/2(J y J z + J z J y ), O zx = V3/2(J Z J X + J X J Z ), 
and O xy = y/3/2(J x J y + J y J x ). Our investigation tells that: First, there is a realistic chance 
to experimentally detect the Thomson scattering signals in TmTe. Second, the intensities 
show different tendency in both materials at several forbidden Bragg spots, which may be 
attributed to difference of the component of the order parameters. 

This paper is organized as follows. Section 2 is spent to introduce a theoretical framework 
in order to calculate the Thomson scattering intensity. In §3, we briefly summarize the CEF 
scheme concluded from our previous paper and explain the ground state both in the absence 
and presence of the applied magnetic field. In §4, we derive some properties of the Thomson 
scattering intensities from the AFQ phase in TmTe, for instance, its dependence on the direc- 
tion of the scattering vector and applied magnetic field. A comparison of the present results 
and those obtained for CeB6 is also found. Finally, §5 is devoted to concluding remarks. Note 
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that a very early stage of the present work is published elsewhere. 17 
2. Scattering Amplitude of Thomson Scattering 



r (6.6')/(G)| 2 , (1) 



The cross section of Thomson scattering is defined as 

In) 

where vq is the classical electron radius. The directions of polarization for the incident and 
scattered photons are denoted by e and e', respectively. The inner product e • e' gives non- 
zero value only when the photon polarization is unrotated; being unity in the a — a' channel 
while cos(2#b) in the it — it' channel where 6-q is the Bragg angle. The scattering amplitude is 
described as /(G) where the scattering vector is defined as G = k' — k with k and k' being 
the wave vectors of the incident and scattered photons, respectively. 

We consider the localized electron system with the (4f) N configuration. The scattering 
amplitude may be given by a sum of the contributions from the localized 4/ electrons: 

/(g) = ^EEE^ 



x (0»(j)\e-^\%(j)), (2) 

where Nq is the number of Tm ion sites. Electron position r n is measured in the coordinate 
system centered at each Tm site j. The |0 M (j)) refers to the ground state of the ./V electrons 
with probability p^(j) where fi distinguishes possible degeneracies. We proceed to rewrite the 
expectation value part in eq. (2), hence we omit the labels j and /x in the following. 

The numerical evaluation of the amplitude can be easily performed by utilizing the so- 
called Rayleigh expansion of the exponential 18 

oo k 

e- G - = 4vr^(-i)^(Gr) £ Y^M^Oa), (3) 

k=0 k z =-k 

where jk means the k-th. order spherical Bessel function and Gr = |G||r|. The solid angles 
of r and G are represented as £1 and Qg, respectively. Note that the similar treatments are 
found in the literatures analyzing Thomson scattering of X-rays from the ordering phase in 
CeB6 for ^-configuration. 16 ' 19,20 

In expanding the ground state, we employ the total angular momentum basis involving 
the radial part, | J, J z ), as follows: 

|0> = £ "WW- ( 4 ) 
■],=-■] 

where we denote \J,J Z ) as \J Z ). The state \J Z ) can be expanded by means of the Slater 
determinant constructed by the one-electron spin orbitals for N electrons. Generally, the 
evaluation of eq. (2) from the iV-electron Slater determinant is tedious, 21 and one can employ 
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the formalism on the basis of the Stevens operator equivalence method. 22 When N = 1 and 13, 
however, the situations are quite simple and we can carry out the evaluation easily. Reflecting 
the fact that Tm 2+ ion is in the /^-configuration, we restrict N = 13 in the following. 

Then, the Slater determinant for thirteen electrons is specified by the quantum numbers, 
orbital (£ z ) and spin (s z ) angular momenta, for a single hole which is the lone unoccupied 
one-electron spin orbital in each determinant. Hence the state \J Z ) is written in the form of 

i s 

\J Z ) = C(J, J z '■ —V-z'i s , — s z )\££ z , ss z ), (5) 

l z =-ls z =-s 

where C(JJ Z : ££ z , ss z ) is the Clebsch-Gordan (CG) coefficient with £ = 3 and s = \ for an / 
electron. The ket \££ z ,ss z ) stands for the Slater determinant for thirteen electrons labeled by 
the hole quantum numbers £ z and s z . Note that when the ket means one-electron spin orbital, 
the minus signs in the CG coefficient of eq. (5) disappear. The one-electron spin orbital is 
described by the product of radial part i? 4 j(r), angular part Y^, and spin part Xs,s z - By 
combining this and eq. (4) with eq. (5), we can continue the evaluation of the expectation 
value of eq. (3), the detail of which is relegated to Appendix. The result is summarized as 

f;(0|e-* G ^|0> = 14<jo(G)> - £ a*(J z )a(J' z )f Jz:J ,. (6) 

n=\ J Z ,J' Z 

Here, we have introduced the amplitude matrix as 

fj„J' z = E C(JJ z :££ z ,ss z )C(JJ' z :U' z ,ss z ) 

&z $z 

OO 

k=0 

x c k (££ z ,££' z ), (7) 
where the Gaunt coefficient is defined by 

= V^TT / y^M^-i'^Yia^ ( 8 ) 

and 

poo 

(j k (G)) = / r 2 j k (Gr)Rl f (r)dr. (9) 
J o 

Since the first term in eq. (6) is independent of the ground state, it has no contribution to 
the scattering intensity as far as G is chosen as antiferro-type spot. 

The Gaunt coefficients are evaluated by means of the Wigner 3j symbols. In the present 
case of / electron system with £ being fixed to three, only the terms for k = 0, 2, 4, and 6 are 
relevant. As a consequence, fj zt j' and eventually the scattering amplitude itself are invariant 
under the transformation G <-> — G. In this context, we do not discriminate between G and 
— G in the present work, in particular, when we perform numerical evaluation of the scattering 
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Fig. 1. (Color online) Radial integrations of spherical Bessel function (jfc(G)) for divalent Tm ion 
as a function of G = |G|. The bold solid (black), bold dotted (black), thin solid (red), and thin 
dotted (red) lines for k = 0,2,4, and 6, respectively. 



intensities. Notice that we can verify the symmetry relations exhibited by the amplitude matrix 
elements as follows 

fj„JL = fj'M = i-) Jz ' Ji f-J z ,-J' z = {-) J *- J *f-J>,-J z . (10) 

A remaining task is to calculate the coefficient a(J 2 )'s. We briefly summarize the methods 
and results in the next section. 

3. Ground State 

Thulium telluride (TmTe) is a magnetic semiconductor crystallized in a cubic, NaCl struc- 
ture with a lattice constant of a = 6.35A. The Tm ion is in a divalent state with one 4/ hole 
[(4/) 13 ] configuration ( 2 F7_). The radial part of the wave function i? 4 j(r) we use is calculated 
within the Hartree-Fock approximation for Tm 2+ . 23 The (jk{G))'s are evaluated by means 
of i?4j(r) as shown in Fig. 1. Then, the angular part of the wave functions are prepared as 
follows. 

Under the cubic CEF potential, the ground multiplet spanned by J = | subspace is split 
into two doublets Tq and and a quartet T$. Since their total separation is believed to be 
around 15 K, 12 we should retain all the bases. In the previous paper, we have introduced a 
model Hamiltonian on the J = 7/2 multiplet basis to describe the phase diagram for TmTe. 15 
Analyzing carefully an interplay among the CEF potential, the Zeeman energy, and multipolar 
interactions, we have concluded that the T3 AFQ order parameter and a CEF level structure 
r8-r6-r7 naturally explain the observed field dependence and anisotropy of the phase diagram. 
In the present analysis, the intensity of X-ray scattering is calculated by using the mean- field 
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ground state derived from the same Hamiltonian. Here, let us summarize the Hamiltonian 
and its mean-field results briefly. 

The basic assumption in the model is that the original fee lattice of Tm ions can be 
decoupled to four distinct sc sublattices. 24 Then, the model on the sc lattice is defined by a 
sum of three parts, H = He + Hz + Hq, where 



Hq is the CEF Hamiltonian defined in ref. 25, and Hz is the Zeeman energy in the magnetic 
field H with g = 8/7 being the Lande g factor for J = 7/2. The T3 AFQ interaction relevant 
for TmTe is given by Hq, where the summation over (ij) is restricted to the nearest-neighbor 
sites in the sc lattice. For simplicity we do not consider influences of field-induced multipoles 
in this model Hamiltonian. 

Concerning the parameters in the CEF Hamiltonian Hq, we assume W = — 0.417K and 
x = 0.5 which lead to a level scheme r8(0K)-F6(5K)-F7(10K). This is nothing but scheme (b) 
in ref. 15, namely the most promising CEF level scheme for TmTe. The quadrupole coupling 
constant is determined so as to give the transition temperature Tq = 4K at zero field, for the 
fixed CEF level scheme. We expect that the transition temperature must be suppressed and 
becomes closer to the real value Tq = 1.8K when the strong fluctuation is taken into account. 

Applying the mean field approximation for the AFQ interaction, one can determine the 
stable order parameters depending on the direction of the magnetic fields. 15 It is shown that 
the 0° order appears in H \ | (001) whereas the 0\ order is stabilized in H \ | (110). On the other 
hand, Of> an d 0\ are almost degenerate at zero field and in H ||(111) due to high symmetry. 
Thereby, we assume the 0\ order at zero field in this study, because the observed field- 
induced antiferromagnetic structure in H ||(110) indicating 0\ are continuously connected to 
the zero field. 13 These mean field analyses obviously provide the ground state wave function 
at each sublattice site, which can be used to calculate the X-ray scattering intensity at zero 
temperature. 

Finally, we shall briefly comment on the properties of the AFQ domains within the model 
described by eqs. (11a) ~ (He). As discussed above, the model leads to a simple antiferro- 
type structure characterized by a wave vector Ki = (111) in units of ir/a. Although this wave 
vector is unique on the sc lattice, it allows degeneracy on the fee lattice, with K2 = (HI), 
K3 = (111), and K4 = (111). This degeneracy is equivalent to the degeneracy when combining 
four decoupled sc sublattices to a single fee lattice. Therefore, the four K domains remain 




(11a) 




(lib) 




(11c) 



(ij) 



6/18 



J. Phys. Soc. Jpn. 



Full Paper 



to be unchanged in the present model even when the magnetic field is applied. The stability 
of the K domains is determined exclusively by a subtle inter-sublattice interaction. 26 In the 
present paper, we will present the results of each K domain and will not discuss details on 
the stability problems. 

4. Thomson Scattering Intensities 

4-1 Remarks on scattering from K domain 

The scattering vector in Thomson scattering to detect the antiferro-type ordering pattern 
in TmTe is simply described as G = (2h + 1, 2k + 1, 2£ + 1) with h, k, and I being integers, 
in units of it /a. In this case, the phase factor appeared in eq. (2) becomes e _lG ' R J = +1 or 
— 1 depending on which sublattice Rj belongs to. We can verify that for any antiferro-type 
G, there exists only one m among 1 ~ 4 which satisfies q~ iG ^j = e -*K m -Rj a ^ ever y j^. 
Thus, the scattering amplitude remains finite only from the if-domain which satisfies this 
relation. In this sense, the scatterings from distinct K-domains should be identified by the 
corresponding scattering vectors. 

Here, we calculate the intensity for the perfect single .fT-domain, which is picked up by 
the scattering vector. This should be kept in mind when we compare the calculated results 
with the experimental ones. If each i^-domain would have nearly the same population, our 
results overestimate factor four. 

4-2 Dependence on the direction of scattering vector 

We exploit a group theoretical analysis on the scattering amplitude /(G) [eq. (2)]. The 
amplitude is invariant under the symmetry operations keeping both crystal and G unchanged. 
Under cubic symmetry (Oh), /(G) is constructed by the quantities belonging to the identical 
(ri) representation in a point group Oh x G. Here, some symmetry operations in Oh are 
forbidden by assuming an artificial strain along G in Oh x G. Since /(G) is expanded by 
a linear combination of the terms with even rank, only the even rank multipole operators 
belonging to the T\ representation inO/,xG contribute to the scattering intensity. This is a 
striking difference compared with a starting point of the similar analysis for the magnetic 
neutron scattering form factor where the unprojected scattering operator behaves as the 
odd rank multipole operators. 27 In the present case of the antiferro-type ordering phase, it 
corresponds to detect the T\ representation from rank two operator. This is easily confirmed 
by re-expressing eq. (6) as 



13 



J>|e- iG - r "|0) = 130o(G)) 



n=l 




(12) 



fc=2,4,6 J Z ,J' Z 
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where coefficient B k (J z , J' z ) is obtained by eliminating expression ^2^Lo(jk(G)) from the right 
hand side of eq. (7). The first term is canceled by the contributions from two sublattices, which 
leaves the expansion starting with the term k = 2. 

Equation (12) indicates the presence of the contributions from the terms of rank four 
and six. However, qualitative behavior of the whole intensity is well understood by that 
from the leading term proportional to (.72(G)). The reasons are two fold. First, because the 
symmetry properties of the coefficients B 2 (J Z , J' z ), B 4 (J Z , J' z ), and B & (J Z ,J Z ) deduced from 
relation similar to eq. (10) are the same one another, the latter two terms do not give rise 
to qualitatively new properties which are absent for B 2 (J Z ,J' Z ) term alone. It means their 
influence on the total intensity is quantitative, not qualitative. Second, it turns out from 
the numerical calculations in the following subsections that the contribution from the term 
proportional to (.72(G)) dominates the intensity. Therefore, though our numerical calculation 
shall include the contributions from the terms with rank four and six, we proceed to make a 
group theoretical consideration deduced only from the rank two term and derive some selection 
rules which qualitatively explain the behavior of the whole intensity. 

In a point group Oh, rank two quadrupole operators are T3-type (O^O^) an d I^-type 
(O yz , O zx ,O xy ). For G || (001), (111), and (110), the components to be invariant under Oh x G 
are O®, (O yz + O zx + O xy )/-\/3, and O® and O xy , respectively. Note that this is the same as 
symmetry lowering of quadrupoles by the magnetic field, as discussed in ref. 28. Although 
antiferro-type G spot does not exist in the (001) nor (110) directions, they are interpreted 
as the limiting cases of the spots, for example, at G = (1, 1, 2h + 1) and (2h + 1, 2h + 1, 1), 
respectively. 

In the absence of the applied field, the order parameter is 0\ as explained in the previous 
section. In the cubic symmetry, two more independent primary order parameters are obtained 
by rotating O 2 by an angle ±27r/3 about the wave vector specifying the K-domain. The 
domains specified by these primary order parameters are called as S-domains. In the K\- 
domain, for instance, the primary order parameters of three S'-domains become 0\, ( K \[?>0\ — 
0^/2, and — (y/SO^ — O^)/^- Then, for a given G, if the Ti representation contains 0° 
and/or 0\, we can expect the Thomson scattering intensity remains finite. We proceed to our 
investigation assuming that three S'-domains have the equal population in each i^-domain. In 
the following, the numerical results are presented for the a — a' channel when unspecified. 

Let us consider G||(lll) at first. It is clear that the T\ representation (O yz +O zx +O xy ) / '\/3 
in Oh x G is not involved in the order parameter. Therefore, the scattering is forbidden in this 
case. Then, we examine the intensity at G = (2h + 3, 2h + 1, 2h + 1), which approaches to G || 
(111) in the limit of h — > 00. As shown in Fig. 2, the intensity obtained at (2h+3, 2h+l, 2h+l) 
is very tiny as expected. 

Second, we consider G = (001). Two of the three S-domains whose order parameters 
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Fig. 2. (Color online) The Thomson scattering intensities per Tm ion under no applied held in the 
a — a 1 channel. The (black) circles, (red) diamonds, and (blue) triangles are the intensities for 
G = (2/i+l, 1, 1), (2/i+l, 2/i+l, 1), and (2/i+3, 2/i+l, 2/i+l), respectively, with h = 0, 1, 2, • • • . The 
(black) solid and (red) dotted lines show the limiting curve for G || (100) and (110), respectively. 

include can give finite intensities in this case. Obviously, the same result is expected for 
G = (100) and (010). In Fig. 2, we plot the calculated intensity at (2h + 1, 1, 1), together with 
the result for continuous G || (100). The latter is evaluated by assuming e ~ lG '^ in eq. (2) 
being +1 or —1 corresponding to R/s sublattice. As seen from Fig. 2, the limiting curve is in 
good accordance with the intensities obtained at the real antiferro-type spots G = (2/i+l, 1, 1) 
even when h is small. 

Then, for another limiting case, G || (110), the discussion similar to that for (100) is easily 
confirmed. That is, the direction (110) is considered as the limiting direction of an antiferro- 
type spot G = (2h + 1, 2h + 1,1). The scattering intensity at (2h + 1, 2h + 1, 1) is equal to 
those at the corresponding spots at (1, 2h + 1, 2h + 1) and (2h + 1, 1, 2h+ 1). We verify these 
results and display the curve together with that for G || (110) also in Fig. 2. 

4-3 Dependence on the direction of applied field 

When external field is applied to the system, usually, degeneracies associated with the S- 
domains are lifted depending on the direction of the field, which also removes the equivalence 
of the intensities under no external field for several high-symmetry G directions. For instance, 
for H || (001), the primary order parameter of the ground state becomes O®- The identity 
representation in O h x G is 0%, (\/30| - 0£)/2, and (-\/30f - 0° 2 )/2 for G || (001), (100), 
and (010), respectively. Then, we expect the intensities for G || (100) and (010) are equivalent 
and weaker than that for G || (001). For corresponding antiferro-type G = (1, 1, 2h + 1), (2h + 
1,1,1), and (1, 2h + 1, 1), the intensities for the latter two are the same, which are weaker than 
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Fig. 3. (Color online) The Thomson scattering intensities per Tm ion under the applied field along 
(a) (001) and (b) (110) directions in the a — a' channel. The (black) filled and open circles are 
the intensities for G = (2h + 1, 1, 1) and (1, 1, 2h + 1), respectively, with h = 0, 1, 2, • • • . While 
the (red) filled and open diamonds, and (blue) filled triangles are those for G = (1, 2h + 1, 2h + 
1), (2h + l,2h+ 1, 1), and (2h + 3,2h + 1, 2h + 1), respectively. The (black) solid and dotted lines 
show the limiting curves for G || (100) and (001), respectively. The (red) solid and dotted lines 
show the limiting curves for G || (011) and (110), respectively. 



that for the former one as displayed in Fig. 3 (a). Next, we consider the limiting directions 
G || (110), (011), and (101). Since the identity representations belonging to the T3 block for 
them are equivalent to those for G || (001), (100), and (010), respectively, the same relations 
hold. That is, the intensities for G || (011) and (101) are the same, which are weaker than 
that for G || (110). Similarly, the intensities at G = (1, 2h + 1,2ft + 1) and (2ft + 1, 1, 2h + 1) 
spots are the same, which are weaker than the one at (2h + 1, 2h + 1, 1) spot. The intensities 
at G = (2h + 3, 2h + 1, 2h + 1) are negligible and we omit them from Fig. 3 (a). 

For H || (110), the primary order parameter of the ground state is 0\. From the analysis 
for the limiting cases, the series G || (100) and (010) include 0\ as the Y\ representation 
while G || (001) does not. Then, the intensities at G = (2h + 1, 1, 1) and (1, 2h + 1, 1) are 
expected to be the same while tiny intensity, if any, is brought about at G = (1, 1,2ft + 1). 
These tendencies are confirmed numerically as shown in Fig. 3 (b). Similarly, the intensities 
at G = (1,2ft + 1,2ft + 1) and (2h + 1, 1,2ft + 1) give the same values while that at G = 
(2ft + 1, 2ft + 1, 1) is essentially zero. 
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Two / electron systems CeBg and TmTe share some apparent similarities. First, both 
materials are cubic systems. Second, they are considered as the one /-particle systems: In 
CeBg, Ce 3+ ion is in the f 1 configuration in the electron picture, while in TmTe, Tm 2+ ion 
is in the f 1 configuration in the hole picture. Third, they both exhibit AFQ ordering phases 
below the critical temperature. On the basis of these nominal resemblances, our main focus 
is to find the differences they may show. One obvious difference is the J value. Owing to the 
Hund's rule, J = § and \ in CeBg and TmTe, respectively. Under the cubic circumstances 
and inferred from the CEF splitting, their ground states are spanned by one Tg quartet in the 
former, and two doublets and Fj) and one Tg quartet in the latter. Those differences may 
reflect on the differences of the Thomson scattering amplitude between TmTe and CeB6- Due 
to the difference of the value of J and corresponding difference of the bases used, CeB6 does 
not include the term proportional to (je(G)) while TmTe does. 

Other than this difference, their differences tend to be quantitative ones. Among them, 
we examine the experimental results presented by Yakhou et al. They reported the ratios 
of the Thomson scattering intensities at G = (511) and (711) to that at G = (533) from 
the AFQ phase in CeB 6 under no applied field. 4 The results are 1^,(511)//^, (533) ~ 
i^ 3 ^/ (711)//^^., (533) ~ 0.02 in the -it — tt' channel. Here, the scattering intensity in the \i — // 
channel at scattering vector G = (2/j+l, 2/c+l, 2£+l) is denoted as I fl - fl i(2h+l, 2k+l, 2£+l). 
With and without the superscript 'exp.' distinguish between the experimental data and the 
theoretical ones, respectively. The intensities measured in the tt — tt' channel involve the factor 
cos 2 (2#b), which depends on the photon energy. When we compare the experimental data for 
CeBg with those obtained from TmTe, it is convenient to eliminate this factor, which leads 
to the ratios expected from the measurement in the a — a 1 channel. Then, the experimental 
data are interpreted as 1^,(511)//^, (533) ~ 1.22 and 7^,(711)//^, (533) ~ 0.01 in the 
a — a' channel. 

To begin with, we comment on a possibility of the experimental detection of the Thom- 
son scattering signals from TmTe. The strongest signal in Yakhou et aVs data for CeBg is 
obtained at (533). 4 Since they did not present the absolute value of /^_ P T ,(533), we interpret 
the theoretical intensity at this spot for CeBg is strong enough to be detected experimentally. 
Then, our previous evaluations for CeBg correspond to (533) = 1.72 x 10~ 3 from O zx 
and O xy phases, while I c - C > (533) = 3.26 x 10~ 3 from O yz phase in the absence of the external 
field. 16 These values are measured in units of r 2 . per Ce ion site. Although we must take the 
population of the S'-domains into account when we compare the theoretical results with the 
experimental ones, it may be reasonable to claim the intensity around 1.0 x 10~ 3 is detectable. 
In the same units, our present numerical calculation tells that I a -a< (533) = 5.63 x 10~ 3 for 
TmTe. Intensities at another spots such as (511) and (711) give three to eight times larger 
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than that at (533). Thus, we assert that the experimental detection of the Thomson scattering 
intensity in TmTe is realistically attainable. 

Next, we investigate the ratios. For CeB$, the magnitudes of the calculated intensities 
are compatible with the tendency in the experiments. That is, /„-_„-/ (533) is several orders 
of magnitude stronger than 7 7r _ 7r /(511) and 7 7r _ 7r /(711). Precisely, I 7r _ 7r /(511)// 7r _ 7r /(533) ~ 
1(T 3 and 7 7r _ 7r ,(711)// 7r _ 7r /(533) ~ 10~ 2 in our calculation 16 while 1^,(511)//^, (533) ~ 
I^ p n i (711) / I^ w , (533) ~ 10~ 2 in the experiment. 4 In a qualitative sense, we believe the dif- 
ference at G = (511) is irrelevant considering the given circumstances such as a lack of 
information on the domain population, the weak signals at G = (511) and (711), and so on. 
Our calculations show i ._ CT /(511)/i CT _ (T '(533) ~ 3.01 and l a - a i (711)// cr _ -' (533) ~ 8.14 from 
the 0\ phase in TmTe. 29 That is, 1*^., (711) is much smaller than 7^^., (511) and -^^(533) 
for CeBg, while these three quantities have nearly the same magnitudes for TmTe. We believe 
this difference is easily recognized if the measurements are available in TmTe. 

Note that the difference may be attributed to that of the nature between the Ts-type 
order parameters and the T3-type order parameters. Actually, we obtain one corroborating 
evidence that the ratios obtained from CeBg assuming one of T3-type order parameters 0\ 
become I a - a /(511)/I a - a i (533) ~ 6.07 and / (7 _ cr '(711)// cr _ -'(533) ~ 5.68, which are similar 
to the TmTe's values. Hence our concern is to understand why I CT _ (T /(711) gives much larger 
value in the T3-type states than that in the Ts-type states. To this aim, we invoke the group 
theoretical consideration developed in §4.2. The series G = (2h+l, 1, 1) approaches to (100) in 
the limit h — > oo. The intensity of the latter is equivalent to that of (001). Since the scattering 
vector G || (001) detects O® as the Ti representation, finite intensities are expected if the 
primary order parameter of the ground state includes the component O® among quadrupole 
operators. As explained in §4.2, the ground states under no applied field really include the 
O® as the primary order parameter if we take into account all the three 5-domains. Because 
J cr _ (T /(711) is very close to the limiting curve for G || (100) as seen from Fig. 2, J cr _ (7 /(711) 
(with h = 3) may already exhibit the property of the limiting curve. On the other hand, 
the primary order parameter of the ground states of CeBg in the absence of the applied field 
consists of the linear combination of the r 5 -type components. Since they do not involve OS>, 
no intensity is expected in the limit of h — > oo. Thus the larger the value of h is, the smaller 
the intensity becomes, which explains why the intensity at G = (711) is extremely small in 
the ground state of Ts-type order parameter. 

There is one remark on the discussion in the previous paragraph. If our justification of why 
la-a' (71 l)/I a - a i (533) is so tiny in CeBg could be correct, we wonder why the same is not true 
for J cr _ -/(511)// cr _ -/ (533) whose experimental value is ~ 1.22. In our calculation, the ratio 
J cr _ (7 /(511)// cr _ -/(533) is one order of magnitude smaller than that reported by the experiment 
as mentioned before. Consequently, our estimate gives I cr _ -'(511)// -_ cr '(533) ~ 0.1, which 
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is consistent with our justification. We cannot find out the reason why the calculated value 
J cr _ cr / (511) / /cr— o-' (533) differs about an order of magnitude from the experimental one in CeBg. 
Since the ratio lT^ a i (511) / I^ a , (533) ~ 1.22 is inferred from the one in the ir — ir' channel, we 
wait for a direct measurement in the a — a' channel, however, this issue is beyond the scope 
of the present work. 

5. Concluding Remarks 

Owing to the extensive efforts to clarify the nature of the AFQ phase expected from TmTe 
below Tq, the knowledge on the magnetic phase diagram has been established. 11-14,24 ' 30 How- 
ever, detailed understandings of ordered phase, such as the component of the order parameter, 
the nature of the microscopic multipolar interactions, and so on, are still rudimentary, which 
should be addressed. As an attempt toward such direction, in this work, we have investigated 
the intensity of Thomson scattering from TmTe in the AFQ phase. We have introduced a 
theoretical framework to investigate the scattering amplitude on the basis of the Rayleigh ex- 
pansion of the exponential part. 16 ' 19,20 ' 22 By taking the group theoretical idea into account, 
we classify the selection rules determined by the orientation of the scattering vector G. 

In the absence of the external field, combining the rules and the domain considera- 
tion, we have obtained some qualitative criteria of the intensity for the orientation of G, 
which determine the absence and/or presence of the intensity, the degeneracy for several 
G orientations, and so on. When we evaluate the actual intensity, however, we need in- 
formation on the ground states of the system. We have utilized the states deduced from 
the J = | multiplet model developed in our previous work. 15 For ground state with 0\ 
being the primary order parameter under no external field, we have checked the three- 
fold degeneracy of the intensities I(2h + 1,1,1) = /(1,2ft, + 1,1) = 7(l,l,2/i + 1) and 
/(2/i+l,2/i+l,l) = 7(1, 2/i+l, 2/i+l) = J(2/i+l,l,2/i+l), while I(2h+1, 2h+l, 2h+l) = 0. 
Note that the degeneracy and the absence of intensity stated here are concluded from the fact 
that the order parameter is 0|, not from the numerical values of the expansion coefficients. 

Then, we have investigated the cases in the presence of the applied field. The field lifts the 
degeneracy on the 5-domain and breaks the cubic symmetry of G orientation. For example, for 
H || (001), the relations 7(2/i+l, 1, 1) = 7(1, 2/i+l, 1) < 7(1, 1, 2/i+l), and 7(1, 2h+l, 2/i+l) = 
I(2h + 1, 1, 2/i + l) < I(2h + 1, 2h + 1, 1) hold. For H || (110), corresponding relations become 
7(2/1+1, 1, 1) = 7(1, 2/i+l, 1) and 7(1, 1, 2/i+l) = 0, and 7(1, 2/i+l, 2/i+l) = 7(2/i+l, 1, 2/i+l) 
and 7(2/i+l, 2h+l, 1) = 0. These relations have been confirmed by the numerical calculations. 

Finally, we have compared our results with those obtained from another AFQ 4/ electron 
system CeBg. 4,16 We have concluded that the Thomson scattering intensities for TmTe at 
several forbidden Bragg spots are experimentally detectable, for instance, at (533), (511), and 
(711) in the a — a' channel. Then, the fact that the magnitudes of the intensities at several 
spots such as (511) and (711) are different significantly between for CeB6 and TmTe may be 
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ascribed to the difference of the components of the primary order parameters between the 
r^-type and ^-type. The discrimination of the component of the order parameter within the 
rV or Ts-types may be achieved by the measurement of the azimuthal angle dependence of 
the RXS intensity. 15 Since our investigation lacks precise numerical information on the weight 
of K- and S'-domains, we should be careful when comparison of our results with the future 
experimental ones will be attempted. 
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Appendix: A derivation of eq. (6) 

In this Appendix, we explain a brief derivation of eq. (6). From eqs. (4) and (5), the 
expectation value of e~ tG ' r " taken by |0) becomes 

f>|e- G -|0> = E E b i,sM,s>Je z ,s,e>, s >, (A-l) 
™ =1 (e z , Sz ) (e' z ,s' z ) 

where 

h z ,s z = ^2a(J z )C(J,J z :£,-£ z ,s,-s z ), (A-2) 

13 

fe z ,s z -.e z ,s' z = ^(« 2 ,^|e- lG ' r "K^4>- (A-3) 

71=1 

We separate eq. (A-l) into the diagonal and off-diagonal parts as follows. 

13 

^(0\e~ lG ^\0)=I d + I od , (A-4) 

71=1 

where 

Id = E \ b ez,s z \ 2 fe z , Sz -.e z , Sz , (A-5) 

(iz,S z ) 

lod = Yl E t 1 ~ 5 ^A 6 s z ,s'M z: s z h' z ,s'Jt z ,s z -£ z ,s' z - 

(e z , Sz ) (e' z ,s' z ) 



(A-6) 



First, we consider the diagonal part. The expectation value between the Slater determi- 
nants is evaluated as, 

13 1 

Id = Yl i 6 ^i 2 Ey3 

(e z ,s z ) n=l 
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E <<MI|e-* G - r "||CO>- (A-7) 

Here the double-bar state ||^", s™)) represents the one electron state of the coordinate r n , not 
the one hole state. Its representation is ((r\£ z , s z }} = R±f(r)Yt^ z (£l)xs,s z - Since the expectation 
value is independent of the electron coordinate after the integration, we can omit it. Then, 
the diagonal part is rewritten as 

!d= Yl \ b iz,s z \ 2 E fi' z ,s' z -i' z ,s' z , (A-8) 
(«,,«,) (Oi)#(<*>«») 

where 

k,s z :e^ = ((^s z \\e- iG - r \\£' z ,s' z )). (A-9) 
Noticing that J2(e z Sz ) \be z ,s z \ 2 = 1> we obtain 

!d = h' z ,s' z :i' z ,s' z - E l^.sj 2 /^:^- (A-10) 

Similarly, after tedious but straightforward calculations, the off-diagonal part is rewritten 

as 

Iod = EE {-f Z ~ iz ~ lb i z ,sM' z ,s'Ji' z ,s' z -.i z ,s z 
(e z ,s z )(£' z ,s' z ) 

+ E \ h ^s z \ 2 k,s z :t z ,s z - (A-ll) 

Eqs. (A-10) and (A-ll) are combined into the following expression. 

E<0|e-* G --|0> = E h**.A 

n=l (<;,«;) 

"EE (A-12) 

(£' z ,s' z ) 

This expression is simplified with the help of several properties of fi Z)Sz :& z ,s' z which we 
derive below. Substituting eqs. (3), (8), and (9) into eq. (A-9), we get 



fc,s M = S a ^V4^^(-i) k y/2k + T{jk(G)) 

k=0 

x Y kA _ e ,(n G )c k (££ z ,££' z ). (A-13) 
For a diagonal part, the summations over l z and s z give 

oo 

E k,s z :£ z ,s z = 2VSQ-i) t ^Tlfe(G)) 

(e z ,s z ) k=o 

e 

x Y kfi {(l G ) E c k (££ z ,££ z ). (A-14) 
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Since Y?t z =-tC k {U z ,U z ) = {21 + l)5 kfi holds, we obtain 

E k,s M = 2(2£ + l)Vto{jo(G))Y 0t0 = 1400(G)), (A-15) 

for £ = 3. Now, we have an expression 

E<0|e- lG - r 10> = 14< J0 (G)> 

n=l 

- E E ( A -!6) 

in place of eq. (A- 12). 

A close look at eq (A- 13) leads to a relation 

ft z ,s z :t' z ,s' z = {-) lz ~ lz ]-l z -s' z -.-l z -s z - (A-17) 

Utilizing this yields the second term in eq. (A- 16) as 

"EE ^-) iz ~ izb t,sM' z ,s' z h' z ,s' z -.t z ,s z 

(e z , Sz ) (e z ,s' z ) 

= ~ E E b -e' z ,-s' z b -tz,-s z ft z ,s z :e' z ,s' z 
(e z , Sz )(e' z ,s' z ) 

= -EE a *( J 2 M4) E E (A-i8) 

x C(JJ Z : U z , ss z )C{JJ' z : U' z , ss' z )~h z ,s z :t> z ,s> z , 

where use has been made of eq. (A-2). Noticing the definition of fj z j> z [eq. (7)], we verify the 
final expression is nothing but the second term of eq. (6). 
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